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We use the optimal fluctuation method to find the density of electron states in- 
side the pseudogap in disordered Peierls chains. The electrons are described by the 
OO ' one- dimensional Dirac Hamiltonian with randomly varying mass (the Fluctuating Gap 

0^ ■ Model). We establish a relation between the disorder average in this model and the 

0^ ' quantum-mechanical average for a certain double-well problem. We show that the op- 

timal disorder fluctuation, which has the form of a soliton-antisoliton pair, corresponds 
^-j. to the instanton trajectory in the double-well problem. We use the instanton method 

' developed for the double-well problem to find the contribution to the density of states 

from disorder realizations close to the optimal fluctuation. 
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I. INTRODUCTION 



> 

ly-^ ■ In studies of disordered systems, one often encounters the problem of finding the density of single- 

00 ' particle states in energy regions where the states can only exist as a consequence of disorder. An example 

. is the density of negatiAfS-energy states of a particle moving in a random potential (the lowest energy of 

' a free particle is zero).trEl Another example is provided by a superconductor with magnetic impunities. 

. The impurities break Cooper pairs and fill the gap with electron states of arbitrarily small energy.a The 

density of such states ("Lifshitz tails") is typically small, as it requires a large disorder fluctuation to 
induce them. In that case the density can be estimated by the probability of the least suppresaecLdisorder 
fluctuation that creates a state of a given energy (the so-called optimal fluctuation method) UU 

In Ref. ^ we applied the optimal fluctuation method to study the density of states and optical ab- 
sorption in disordered Peierls insulators. In the absence of disorder, the electron excitation spectrum 
''O ' in these quasi-one-dimensional materials has a gap which results from the electron-lattice instability 

Ch , and is accompanied by a periodic lattice distortion with the wave vector Q — 2kF^ We considered 

O ■ half-filled chains, in which case the Peierls instability results in an alternation of long and short bonds 

(dimerization). Disorder in the electron hopping amplitudes induces electron states inside the Peierls 
gap, transforming it into a pseudogap. In Ref. ^we showed that the optimal fluctuation for a disordered 
dimerized chain has the form of a soliton-antisoliton pair. It induces two electron states which lie close 
to the center of the pseudogap and which are symmetric and antisymmetric superpositions of states 
' localized near the soliton and antisoliton. Knowing the typical form of the disorder-induced electron 

states we were able to find the absorption coefficient at small energies. These results were also used to 
explain the coexistence of the dimerization and antiferromagnetism, observed recently in the spin-Peierls 
material CuGeOsE 

The probability of the optimal configuration gives, however, only an estimate of the density of states. 
A more careful calculation requires finding the contribution of disorder realizations close to the optimal 
one. In this paper, we perform such a calculation using the close relation between the optimal fluctuation 
method and the semiclassical approximation in quantum mechanics and fleld theory. In the path integral 
version of this annioximation, the most important paths lie close to the "saddle-points" of the action, 
called instantons.Htl In particular, the instantons for the double- well problem are the classical trajectories 
describing the imaginary time motion between the two wells, which corresponds to tunneling and results 
in level splitting. We show how the disorder average of the density of states can be written in the form of 
a functional integral that represents the sum over all paths of a particle in a corresponding double-well 
potential. Then the optimal disorder conflguration becomes precisely the instanton trajectory and the 
calculation of the density of states can be done using standard instanton methods. 

The outline of this paper is as follows. In Sec. || we introduce the Fluctuating Gap Model describing 
disordered Peierls systems and briefly discuss the previously obtained analytical results for the disorder- 



averaged density of states in that model. Next, in Sec. ^we establish the relation between the disorder 
average in the Fluctuating Gap Model and the quantum-mechanical average for a certain double-well 
problem. We find the instanton trajectory corresponding to the optimal fluctuation and perform the 
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"saddle-point" integration for the average density of states. We discuss our results and conclude in 



Sec. |rV|. Two technical points are addressed in appendices. In Appendix A we discuss how the symmetries 
of the FGM Hamilitonian are crucial for the validity of the optimal fluctuation method. In Appendix B 
we discuss the instanton calculation for a different choice of boundary conditions. 



II. THE FLUCTUATING GAP MODEL 

The Fluctuating Gap Model (ECM) is a continuum model describing electrons in disordered quasi- 
one-dimensional semiconductors." IrH In this paper we use it to describe half-filled Peierls chains with 
disorder in the hopping amplitudes of the electrons between the chain's lattice sites. These systems are 
characterized by the Peierls order parameter, A (a;), which is the (continuum version of the) alternating 
part of the hopping amplitudes along the chain, i.e., the difference between the hopping amplitudes on 
nearest odd and even bonds. In the absence of disorder this order parameter U^sp called "dimerization" ) 
is constant along the chain and equals the gap in the single-electron spectrum.B'E3 Disorder in the electron 
hopping amplitudes results in random variations of the order parameter along the chain (which explains 
the name of the model): 

A{x) ^ Ao + vix) . (1) 

Here, Ao is the average value of the order parmeter and r]{x) is the fluctuating part with a Gaussian 
correlator, 

{r]{x)7j{y)) = AS{x ~ y) . (2) 
The single-electron states close to the Fermi energy ep — 0, are described by the wave function 

^(x) - ( ^^^^^ 

where the two amplitudes i)i{x) and ip2{x) correspond to electrons moving, respectively, to the right 
and to the left with the Fermi velocity vp- Since the density of single-electron states does not depend 
on spin projection, we omit the spin index of the electron wave function. The wave functions satisfy the 
one-dimensional Dirac equation, 

H=l^^'^-^ + cjiA{x)^il;{x)=ei}{x) , (3) 

where ui and (73 the Pauli matrices. For the applicability of the continuum description of electrons, the 
random variations in the hopping amplitudes hpje to be relatively small. 

This model was flrst introduced by KeldyshJid It was also considered in the context of the thermo- 
dynamical properties of quasi-one-dimensional charge-tran|S£ip.salts,ll3 and has been applied to study 
the effect of disorder on the Peierls transition temperatiMjEJO as well as the effect of quantum lattice 
fluctuations on the optical spectrum of Peiecls materials. tJilS Recently it was used to describe the phase 
diagram of disordered spin-Peierls systems .13 

In the remainder of this section we discuss the dependence of the spectrum of the single-electron 
eigenstates on the disorder strength. By p{e) we denote the disorder averaged density of states per unit 
length: 

p(e) = i(Tr[5(/i -£)]), (4) 

where L is the chain length. Due to the charge conjugation symmetry of the Dirac Hamiltonian, the 
density of states is a symmetric function of energy: 

p{-e) = p{e) . (5) 

We also introduce iV"(e), equal to the average number of electronic states (per unit length) with energy 
between and some e > 0: 
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N{e) = / de'p{e') 



(6) 



An analytiqaJjfjejcpression for iV(£) was found by Ovchinnikov and ErikhmanHl by means of the "phase 
formaUsm" 



Nis) 



2A 



(7) 



Here, z = vpe/A and v = upAo/A are dimensionless variables, and J^{z) and N^[z) are, respectively, 
Bessel and Neumann functions. 
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FIG. 1. Disorder-averaged density of states p{£) in the FGM for three values of the dimensionless disorder strength: 
g = 0, i.e., no disorder (curve a), g = 0.25 (curve b), g = 4 (curve c). The free electron density of states was set to unity. 

In Fig.0 we plot the average density of states ^(e) for several values of the dimensionless disorder 
strength g = 1/v = A/{vf^q)- In the absence of disorder (A(a;) — Aq), the electron spectrum has a 
gap between the energies e = — Ao and e = +Ao. Disorder gives rise to the appearance of electron 
states inside the gap. The energy dejjxepdence of the average density of states close to the Fermi energy 
(|e| <C Aq) is given approximately bj; 



Pie) 



(8) 



(here and below we assume e to be positive, which is sufficient in view of Eq.(||)). For g < 2 the density 
of states has a pseudogap (the Peierls gap filled with disorder-induced states). For g > 2 the pseudogaa 
disappears and the density of states becomes divergent at e = 0. This is a Dyson-type singularityjHj 
which occurs in the band center for random Hamiltonians with charge conjugation symmetry. 

In Ref . I we showed that the energy-dependence of the average density of states at small energy and 
weak disorder is mainly determined by the weight of the optimal disorder fluctuation f]{x), which is 
defined as the most probable among the fluctuations that induce a state at given energy e. The form of 
the optimal fluctuation is 



r){x) 



-vpK 



i?, 



tanh K{x — xq -\ ) ) — tanh K{x ~ xq ) 



(9) 



where xq and R describe, respectively, the position and the spatial extent of the disorder fluctuation, 
and K is determined by 



vpK ^ Aotanh(i^i?) 



(10) 



Thus, A(x) = Aq + rj^x), plotted in Fig. ^, has precisely the form of a soliton-antisoliton pair discussed 
in Refs. 21 2^ in the context of polaron excitations in the conjugated polymer trans-polyacetylene. The 
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spectrum of electron states for this A (a;) (also plotted in Fig. |^) consists of a valence band (with highest 
energy — Aq), a conduction band (with lowest energy +Ao), and two localized intragap states tp±{x) 
with energies iso{R), where 



£o(i?) 



cosh(i^i?) 



(11) 



The soliton-antisoliton separation R is fixed by the condition eo{R) = e. The two intragap states are 
the symmetcic and antisymmetric superpositions of the midgap states localized near the soliton and 
antisoliton.Efl The energy splitting 2e decreases exponentially with the soliton-antisoliton separation, so 
that for e ^ Aq, 



i? « In 



2Ao 



(12) 



where = ^f/Aq is the correlation length. 









A,; Y 






-R/2\ 


j+R/2 




£p=0 




FIG. 2. The left part of the picture shows A(x) = Aq + fj{x) for the instanton disorder fluctuation (thick line) and 
the electron density |i/>+(a;)p = for the corresponding intragap states (dotted line); the right part of the picture 

shows the spectrum of electron states for the optimal disorder fluctuation. 



For white-noise disorder with the correlator Eq. 
ri{x) is given by 



the weight p [77(2:)] of the disorder configuration 



p[r]{x)] = exp Jc 



dxrj (x) 

Using Eq. (y_2|) we find that the weight of the optimal configuration, 

p [f]{x)] (X £3 , 



(13) 



(14) 



cf. 



gives, indeed, a good estimate for the shape of the density of states inside the pseudogap for 5 ^ 1 
Eq.(|)]. 

We would like to conclude this section by noting that, whileJii. the random potential problem the 
optimal fluctuation induces one state with large negative energyJiTH in the FGM the soliton-antisoliton 
fluctuation Eq.(^) induces two states with energies ±£. The small difference between the energies of the 
symmetric and antisymmetric states is potentially dangerous for the application of the optimal fluctuation 
method to the FGM. The problem arises in the calculation of the contribution of the disorder realizations 
close to the optimal fluctuation: 

r]{x) — fj{x) + 5ri{x) . 

The perturbation 5h = aidr]{x) can, in principle, strongly mix the symmetric and antisymmetric intragap 
states, because of the small energy denominator 2£ appearing in the perturbation series. Such mixing 
would affect the value of the energy splitting between the two states. This problem arises because of 
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the charge conjugation symmetry of the Dirac Hamiltonian Eq.(^, which imphes the symmetry of the 
spectrum of its eigenvalues around e — 0. We address this question in Appendix A, where we show that, 
just because of the symmetry of the Hamiltonian Eq.(|^), the mixing of the symmetric and antisymmetric 
states is small and thus, despite the small energy splitting, the optimal fluctuation method is applicable. 
The result of Appendix A also justifles our calculation of the optical absorption coefhcient presented in 
Ref. |. 



III. FUNCTIONAL INTEGRATION 

The weight of the optimal fluctuation, given by Eq.(|l^, is the main factor that describes the suppres- 
sion of the density of states at small energy. To obtain a full expression for the density, however, one has 
to take into account the contribution of the disorder realizations close to the optimal configuration. In 
this section we perform this calculation using the correspondence between the averaging over disorder 
realizations r]{x) and the quantum-mechanical averaging over the ground state for a certain double-well 
potential. Then the main suppression factor Eq.(|l^), as well as the correction to it, can be easily found 
using well-known methods developed for the double- well problem. 



A. The corresponding quantum-mechanical problem 

The density of electron states averaged over the white-noise disorder can be written in the form of the 
functional integral: 



Drj{x) exp 



Tr 



5{h - e) 



with 



1 /"^/^ 



dxrfix) 



(15) 



(16) 



-L/2 



and h as defined in Eq.(||). The relation between the functional integral and a particular quantum 
mechanical problem will become apparent after we make a substitution of variables that allows us to 
express the disorder fluctuation ri{x) in terms of the eigenfunction '4>{x) of Eq.(|^) with energy e. Due to 
time-reversal symmetry of the Dirac Hamiltonian its eigenfunctions ^{x) can be chosen to satisfy 



(Ti'0(a;) = ijjlx) , 



(17) 



and, therefore, can be written in the form: 



The two real variables w and (/> introduced by 

u = w (cosh(0) -|- zsinh(0)) 

then satisfy. 



(18) 



(19) 



and 



- vp-r- + e cosh{24>(x)) - Aq = t](x) , 
ax 



dhiw(x) . , ,„ , , 

VF = -esinh(20 a;)) 

dx 



(20) 



(21) 



We now use Eq.(20) to express the disorder configuration r][x) in terms of (t>{x). Then S, defined by 
Eq.([l6|), can be written as 
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S[cl>{x)]^SEmx)]+SB , (22) 



where 



and 



SEmx)] '^dx [ «F g ) + (Ao - £cosh(2(/)))^ ) , (23) 



Sb = \vf (2Ao(/. - esmh(20))|+^/^ . (24) 

»S'_B[0(a^)] has the form of the EucHdean actfon of a particle with mass v\ and coordinate </) moving in 
the imaginary time x in the inverted double-well potential, 

t/(^)--^fl-^^V, (25) 



2 V cosh(2v3) 
plotted in Fig. ||. The angle <p in Eq.(^) is defined by 



cosh(2(p) = — . (26) 



m) 



-<p 


+<p 











FIG. 3. The inverted double-well potential (7(0). The instanton trajectory starts at the left maximum —ip and ends 
at the right maximum 

Finally, in the functional integral Eg. (|l5|) , the disorder strength A plays the role of U. Thus, we 
have identified the quantum double-well problem that correanonds to our disorder average. A similar 
correspondence exists also for the random potential problemHj 



B. The instanton solution 



The instanton 0(a;) is the solution of the classical equation of motion in the inverted double-well 
potential with the boundary conditions 

0(±oo) = ±^ . (27) 

This corresponds to infinitely long motion from the left maximum to the right one. The energy on this 
classical trajectory, 

equals zero, which gives. 
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VF-£ = ±{Ao-ecosh{2(f>{x))) . (29) 

As the solution with the minus sign corresponds to the trivial disorder configuration rj{x) — [cf. Eq.(|2 
we choose the plus sign. Equation (p9[) is easily solved, yielding 



where K and R are defined by 



and 



cosh (a; ^ a^o ^ -f )) 



K ^— tanh(2(/j) , (31) 

VF 



KR = 2ip . (32) 



Using Eq.(g^) we find that the instanton Eq.(^) is exactly the optimal disorder configuration given 
by Eq.(p|) and that Eqs. (p6|), (^l|), and ( |3^ ) defining the parameters of the instanton solution coincide 
with Eqs. ( p^ and (pi]). Finally, for the instanton trajectory, the value of S defined by Eq.(p^) is: 

| = ^ = |(2^-tanh(2^)). (33) 



C. Boundary conditions 

We note that the boundary conditions Eq.(^7|) for the instanton trajectory are not, in general, the 
boundary conditions that one usually imposes on the eigenfunction ^p{x). In particular, ip defined by 
Eq.(p6|), depends on the energy e, while, for instance, the boundary -conditions for a finite Peierls chain 
with an even number of atoms in the continuum limit take the formcJ 

Ini(u(±i/2)) , , 

= ±1 (34) 

Re(u(±L/2)) ' ^ ' 

which is equivalent to 

0i?,L = 0|±L/2 = ±oo . (35) 

We discuss the instanton for these boundary conditions in Appendix B. It is clear, however, that in the 
limit L — > oo the density of states should not depend on the boundary conditions. Therefore, we choose 
here the ones that lead to the easiest calculation of the density of states, namely. 



±f ) , (36) 



where 0(x) is the instanton solution given by Eq. (pOj) with .tq = 0. With the last condition <j){x) is 
antisymmetric, so that 

(t>L = -^R . (37) 



D. Regularization 

In order to perform the functional integration in Eq. ([l5|) , we first have to regularize it. We do this by 
dividing the range of variation of x, [— -j, -j], into N small intervals of size a, assuming the value of rj to 
be constant inside each interval. Then, 



7 



1 /-i/a ^ 



and the integration measure, 



Dvix) = n ' (39) 

is defined in such a way that the average of the unity operator is one. 

In the regularized version of the substitution Eq.(|2^), the new variables (/>„ (n = 1, . . . ,N) are the 
values of (j) on the right end of each interval and 

(l)o = (t>L, (40) 

is fixed by the boundary condition at the chain's left end. For the moment we treat Eq.(^0|) as a 
functional substitution and forget about the boundary condition at the right end. Up to second order 
in a, the regularized form of the substitution reads: 

{4'n~4>n-l) (£COsh(2(/)„_i) - Ao - ?7„) , , , 

Vp = T r . (41) 



(l-f^sinh(20„_i)) 



It is necessary to retain the higher order term in powers of a for the correct calculation of the Jacobian 
of the substitution. The integration measure can now be written as 

Dvi.) - n ' (42) 

„=i [27raA] = 



with 



J = TT fl - — sinh(2(^„_i)~^ w exp [ — — [ sinh(2(?;.(a;)) | 



(43) 



It is important that Eq.(^l|) provides a one-to-one correspondence beween the disorder configuration 
{r]n} and the values of the solution of the differential equation ( ^o|) at the points x = an , n = 1, . . . , N . 
In other words, Eq.(^l|) relates the disorder configuration to the electron wave function in the presence 
of this disorder, which gives us a convenient representation for the trace of the delta-dunction in Eq. (|l5| ) . 
The energy e is the eigenvalue of the Hamiltonian ft, if at the right end satisfies the boundary condition: 

(f>N ^(l)R ■ (44) 
Therefore, the average density of states Eq.(p^ can be written in the form: 



p{e) = ^ jDT]{x)e^p ^-^ 



de 



S i^N - 0fl) . (45) 

00, '/(a;) 



The delta-function removes the integration over (p^ and, using Eq.(pO|), one can easily calculate the 
derivative of (j)N with respect to e, keeping both the disorder r](x) and the value of (j) at the left end of 
the chain fixed: 

i^N\ _y(f) /■^/'^^cosh(20(a;)) 



dx \7 " . (46) 



Here, y{x) is given by 



9^ /0o,r,(x) J-L/2 y{x) 



?/(a;) = exp f — / dT'sinh(2(/i(x')) | , (47) 

\VF J-L/2 I 



8 



so that J, defined by Eg. (|4^) , can be written as 

J=^^. (48) 

We thus arrive at the following final form of the functional integral 

p{e) = \ Jo'cPix) exp (^-^) F[<j>{x)] , (49) 
where the integration measure is defined by 

D<j>{x)^ r 11 T' (50) 

[2TraA]^ [2TraA]^ 

and 



E. Saddle-point integration 



In this section we perform the functional integration in Eq.(49) over configurations </>(a:) in the vicinity 
of the saddle-point configuration (l)sp{x) defined as 



SS 



S(l}{x) 







(52) 



<^sp(a;) 



The boundary conditions Eq.(p6|) were chosen to ensure that (j)sp{x) — (t>{x), where (t>{x) is given by 
Eq.(|ol) with xq = 0. 

The result of the integration is: 



Pie) 



JD /d(j)N 

~T [~dr 



exp 



So 



0o.»)(2;) 

where Sq is given by Eq.(p3) and D is the determinant of the operator 



4^ 

2 dx^ 



+ Ae (e cosh(40(a;)) - Aq cosh(20(x))) , 



(53) 



(54) 



obtained by the second variation of the Euclidean action Eq.(g3|). In Eq.(p3D J, D, and d(f>N/de have to 
be calculated at (p{x) = 0(a;). 

Using Eq.(p9|) (with the plus sign) one easily finds that 



— sinh(2,^) = --^ln( ) , 
vf ax > -J- ' 



dx 



so that Eq.(^7|) can be written in the form 

y{x) = 



dx 



(x) 



(55) 



(56) 



Since (j){x) is antisymmetric, y (-j) = 1, so that J = 1 and 



90. 



'JV 



de 



_ sinh (20(f)) 



(57) 
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As for large L, sinh (2(/>(-j)) can be replaced by sinh(2i^) and 



dx 

we obtain 



«2A'smh(2¥.)e-^^, (58) 



On 



(59) 



We calculate the determinant D using the standard method for calculation of determinants of 
Schrodinger operators in one dimension, which can be found in Ref. |^. The result is 

D - L^ll^^ sinh(2^) e-^^ . (60) 

The determinant is proportional to the chain size, because the soliton-antisoliton pair caujappear at any 
place in the chain (formally the factor L comes from the integration over the zero modetl). Taking all 
factors together, we finally obtain for the averaged density of states 

p(e) = ^sinh(2^)exp('-^') . (61) 



nA \ A ^ 

The resulting expression for the average density of electron states in the limit e <C Ao is 

Pie) = (:^) ' ' . (62) 

For 5 <C 1, this agrees with Eq.(||), confirming the validity of the optimal fluctuation method at small 
energies and weak disorder. 



IV. DISCUSSION 



In this paper, we have used functional integration to establish a relation between the Fluctuating Gap 
Model, which describes one-dimensional electron motion in the presence of both a Peierls distortion and 
quenched disorder, and the quantum motion of a particle in a double-well potential. Averaging over 
disorder realizations in the model with quenched disorder corresponds to the sum over all paths in the 
quantum-mechanical problem. We have shown that the instanton trajectory describing the imaginary- 
time motion between the two wells corresponds to the optimal disorder fluctuation. The probability of 
this fluctuation determines the asymptotic behavior of the averaged electron density of states. 

We showed that the most probable form of the wave function of the electron states lying deep within 
the pseudogap contains two peaks. The optimal disorder fluctuation that induces such a state, has the 
form of a soliton-antisoliton pair and the peaks of the wave function are localized near the two kinks of 
this fluctuation (see Fig. ^). As we demonstrated in this paper and in Ref. ||, the instanton approach 
allows for a relatively easy calculation of the density of states and absorption coefficient. Our result 
Eq.(|62[) is valid if the density of disorder-induced states is small, which is the case when |e| <^ Aq and 
5< 1- 

Above, we emphasized the strong relation between the calculation of the disorder average and the 
double-well problem. However, there is also a clear difference: In the double-well problem, an important 
role is plajjed by multi-instanton configurations, in which instantons are followed by anti-instantons and 
vice versus Although the action for the configuration containing n instantons and anti-instantons is 
n times the action of a single instanton 5i„st, the multi-instanton contribution obtains a factor T'^/n\ 
from the sum over all possible instanton positions (T is the imaginary time corresponding to L). This 
statistical factor grows with T, while the suppression factor exp(— nS'inst/?j) does not decrease any further, 
making the multi-instanton configurations important at large enough T. We note that it is the sum over 
all possible numbers of instantons and anti-instantons that gives the energy splitting between the two 
lowest states in the double-well potential.cl 
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By contrast, the asymptotic expression Eq.(62) for the density of states in the FGM is given by only 
one instanton, which means that the largest contribution to the averaged density of states comes from a 
single disorder fluctuation. The origin of the difference lies in the fact that, while in the quantum double- 
well problem the anti-instanton describes a real physical process (the tunneling between the wells in the 
direction opposite to the one described by the instanton), in the FGM it corresponds to zero disorder. 
To see this, wc note that the instanton 0a(2:) — 4){—x), obtained from the instanton by "time-reversal". 



X — > —X, is the solution of Eq.(p9|) with the minus sign. Comparing this equation with Eq.(20), we find 
that 4'a.{x) corresponds to r]{x) — 0. The Euclidean action Se of the anti-instanton is exactly cancelled 
by the "boundary" term Sb [cf. Eq.(p^)]. It is then clear, that anti-instantons (and, hence, the multi- 
instanton configurations) cannot play any role for the calculation of the density of states. Formally, this 
happens because for (l){x) = ipa.ix) we have: 



) vfKs 

which tends to a constant value in the large L limit, while the corresponding expression for the instanton 
solution Eq.(|59|) diverges as exp(ifL). Since the value of the determinant D oc iexp(— ifi) for the anti- 
instanton is the same as for the instanton, the anti-instanton contribution is proportional to L exp(—KL) 
and, therefore, vanishes in the L oo limit. 
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APPENDIX A: THE ROLE OF SYMMETRIES OF THE FGM HAMILTONIAN 

In this appendix we show that owing to the symmetry properties of the FGM Hamiltonian a small 
perturbation of the optimal fluctuation 

r]{x) = f]{x) + Sr]{x) , (Al) 

does not strongly mix the symmetric and antisymmetric states, despite the small energy splitting between 
them. As was explained in Sec. |ll| such a mixing, if it would occur, would invalidate our optimal 
fluctuation calculation. 

First, we obtain the effective Hamiltonian, acting on the subspace of the two intragap states, which 
includes the virtual excitations to all other electron state s. T o this end we write the single-electron 
Hamiltonian for the disorder realization ri(x), given by Eq. (|Al|) , in the form: 

/i = /i(o) + , (A2) 

where 

= ^3 ^ ^ + ai ( Ao + fjix)) , (A3) 

and /i'^-* is the perturbation: 

= ai6rj{x) . (A4) 

We divide the Hilbert space of the Hamiltonian into two subspaces: A and B. The subspace A 
consists of the symmetric and antisymmetric intragap states, which here are denoted as |±), while all 
the other states belong to the subspace B, so that 

Pa + Pb^I, (A5) 
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where Pa,b are the projection operators for the corresponding subspaces. The effective energy-dependent 
Hamiltonian, acting on the subspace A, has the form: 

h^{e) ^ h^^^ + h'^^^ + h'-2\e) (A6) 
where the first term is the diagonal 2-by-2 matrix: 

p^/.(")p^ = f +^°(^) (A7) 



-eo(i?) 
the second term is the first-order perturbation, 

h^^^ = pAh^^^pA , (A8) 



and the last term, 



describes the higher-order processes involving virtual excitations to the subspace B. 

We now want to show that the non-diagonal matrix elements of the effective Hamiltonian ^.^(e) 
are smaller than its diagonal matrix elements. First, we show that the first-order nondiagonal matrix 
elements equal 0: On the one hand, using the relation '4'±{x) = <7i'ip±{x) [cf. Eq.(p^)] we obtain that the 
first-order perturbation Hamiltonian is real and symmetric: 

(-f|/iW|-) = = {+\h^'^\-r = ■ (AlO) 

On the other hand, we can use the relation between the positive-energy and negative-energy eigenstates 
of the Dirac Hamiltonian, ip±ix) = ^o'2V'^(a^), to show that 

i+fh^'^n = -{+W2h^'^<^2\-) = -{-\h^'^\+) ■ (All) 
Comparing Eq.(AlO) with Eq.( All[ ) we obtain: 



= =0. (A12) 

Next we consider the non-diagonal matrix elements of h'"J^\e). The perturbative expansion of h'"j\e) 
in powers of /i^^-* describes virtual excitations into the subspace B, so that all intermediate states in this 
expansion are separated from the intragap states by a large energy eo{R)). Since |e| ~ £o{R), the 
expansion terms involved do not contain small energy denominators. Moreover, the non-diagonal matrix 
eleme nts of the operator ( [A9| ) equal at e = 0. The proof of this statement is analogous to the proof 
of Eq.(A12) and we shall not repeat it here.|^ Hence, the non-diagonal matrix elements of the operator 



h^2\e) are 0(e). Moreover, they contain an additional small factor Sr]{x) /Aq. Therefore, the small 
disorder perturbation 6ri(x) cannot result in a strong mixing of the two intragap states. 

Next we show that the corrections to the energies of these two states due to the diagonal matrix 
elements of the effective Hamiltonian are also small, which will justify the validity of our calculation of 
the optical absorption coefficient .0 Using the expressions for the wave functions of the symmetric and 
antisymmetric statesjj the straightforward calculation of the first-order terms gives 

(+|/,(i)|+) = = JdxSv{x)cb{x - ^)cf>{x + |) , (A13) 

where 



^The important point is that in the basis of the eigenstates of ft'"', satisfying Eq.(^), ft'"' and ft'^-* are real 
symmetric matrices and that they both anticommute with a2- 
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Since the condition eo{R) ^ Aq implies R ^ [see Eq.(^], the last equation can be written in the 
form: 

« le-t r\x5^{x)^e,{R) , (A15) 

so J-B./2 J-R/2 ?0 ^0 

i.e., the first-order correction to the energy of the intragap state with the energy ±£o(^) is proportional 

to this small energy times a factor O (^^^^^ ■ This result can be easily traced back to the fact that the 

energy of the midgap state of a single kink cannot be changed by a small perturbation of A(x) (see, e.g., 
Ref. pq) and, therefore, it holds in all orders of the expansion of the effective Hamiltonian in powers of 

Thus, we found that, despite the small energy splitting between the symmetric and antisymmetric 
intragap states, a small variation of the disorder fluctuation does not strongly affect either the wave 
functions, or the energies of these two states. This result is a direct consequence of symmetries of the 
Hamiltonian Eq.(0). 



APPENDIX B: DIFFERENT BOUNDARY CONDITIONS 

In this appendix we discuss the instanton trajectory for the boundary conditions Eq. (p5|) , which one 
obtains in the continuum limit for an open chain with an even number of atoms. To avoid spurious 
divergences related to an infinite value of (f> at the boundaries, we will use the following procedure. First, 



we impose, instead of Eq.(35), the boundary conditions, 

<f>R:L = 0l±L/2 - ' (Bl) 

where $ 3> V' is large but finite. Next we take the limit L oo and only then (f> is tended to infinit y. 

Since for |$| > part of the instanton trajectory 0(a;)lies outside the well of the potential Eq.(p5|) 
located between —(f and (/?, the energy E defined by Eq.(p8) is positive. For E > the equation for the 
instanton trajectory is. 



= y^(Ao - £ cosh(2<?!)(a;)))^ + 2E . (B2) 

Now, the particle's velocity becomes very large as soon as it leaves the well. The large Euclidean action 
Eq.(p3|) associated with this region is, however, completely compensated by the surface term Eq.(p^ 
in the infinite L limit. To see this, one notes that even for infinite $ the time of motion L along the 
instanton trajectory is finite, as the velocity grows very fast outside the well. We are, however, interested 
in long time (large L) trajectories. The only way to make L large is to tend the positive energy E to 
zero, in which case the particle would spend most of the time near the maxima —tp and ip of the potential 
U{(j)), where it has a small velocity. The relation between the energy and the time of motion for large L 
is, 

E^^^-^^e-^^e-^"^ . (B3) 

For ^ If the dependence of the energy on <i> can be neglected. 

When the energy E tends to zero, the trajectory given by Eq.(B2) inside the well, i.e., for \(j)\ < ip, 
tends to the solution of Eq.(p9|) with the plus sign. In other words, in the limit L — > oo, this part of 
the trajectory becomes the instanton solution Eq.(|30|) that we considered above. On the other hand. 



outside the well (|0| > <p) the solution of Eq.(B2) tends to the solution of Eq. (|29|) with the minus sign. 
The latter corresponds to zero disorder, r\(x) — 0, which explains the cancellation between the Euclidean 
action coming from outside the well and the surface term. We proved, the refo re, that in the limit L 



the suppression factor S for the instanton with boundary conditions Eq.(Bl) is identical to Eq.(p3| 



The boundary conditions Eq.(Bl) imply that the solution (\){x) is antisymmetric. Therefore, the 



Jacobian J given by Eq.(}48|) equals unity, as before. At the same time, however, the values of the partial 
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derivative ( and the determinant D differ considerably from Eqs. ( ^9|) and (|60|), respectively. 



For instance, 



5e ^^_^(,) 16sinh4(2¥.) ' ^^^^ 

Nevertlieless, the product of the determinant and the partial derivative remains unchanged up to small 
corrections that vanish for L — > oo. Therefore, despite the fact that the instanton trajectories for the 



boundary conditions Eq.(|35|) and Eg. (|36D are rather different, the average density of states obtained by 
the saddle-point integration near these trajectories is the same for sufficiently long chains, as, of course, 
should be the case. 
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